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Abstract. In this paper we present a modified version of the CORS method based on a new calibration of the Surface Brightness 
function in the Stromgren photometric system. The method has been tested by means of synthetic light and radial velocity 
curves derived from nonlinear pulsation models. Detailed simulations have been performed to take into account the quality of 
real observed curves as well as possible shifts between photometric and radial velocity data. The method has been then applied 
to a sample of Galactic Cepheids with Stromgren photometry and radial velocity data to derive the radii and a new PR relation. 
As a result we find logR = (I.I9 + 0.09) + (0.74 ± 0.1 1) log P (r.m.s=0.07). The comparison between our result and previous 
estimates in the literature is satisfactory. Better results are expected from the adoption of improved model atmosphere grids. 
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1. Introduction 

Classical Cepheids are the cornerstone of the extragalactic dis- 
tance scale. Thanks to their characteristic Period-Luminosity 
(PL) and Period-Luminosity-Color (PLC) relations they are tra- 
ditionally used to derive the distances to Local Group galaxies, 
and (with the advent of space observations) to external galax- 
ies distant up to about 25 Mpc (targets of a Hubble Space 
Telescope Key Project, see Freedman et al. 1997, 2001). As 
primary indicators they are used to calibrate a number of sec- 
ondary distance indicators (see e.g. Freedman et al. 2001) 
reaching the region of the so called Hubble flow where the 
Hubble law can be applied and an estimate of the Hubble con- 
stant can be derived. 

Moreover, the comparison between Cepheid physical parame- 
ters (stellar mass, luminosity, chemical composition) based on 
evolutionary and pulsation models supplies the unique opportu- 
nity to pin point the occurrence of deceptive systematic errors 
(Bono et al. 2001a; Moskalik 2000) on the Cepheid distance 
scale. 

In particular radius determinations are important to constrain 
both the intrinsic luminosity, through the application of the 
Stefan-Boltzmann law, provided that an effective temperature 
calibration is available, and the stellar mass, by adopting a 
Period-Mass-Radius relation (e.g. Bono et al. 2001). 
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Many investigations have been devoted during the last decade 
to the derivation of accurate Period-Radius (PR) relations for 
Classical Cepheids both from the empirical (see e.g. Laney 
& Stobie 1995; Gieren, Fouque, & Gomez 1998; Ripepi et 
al. 1997) and the theoretical (Bono, Caputo, Marconi 1998; 
Marconi et al. 2003) point of view. 

Empirical Cepheid radii are generally derived either by means 
of the Baade Wesselink (BW) method (Moffet & Barnes 1987, 
Ripepi et al. 1997; Gieren et al. 1998, just to list a few exam- 
ples) both in the classical form and in subsequent modified ver- 
sions, or with interferometric coupled with trigonometric paral- 
laxes techniques (Nordgren, Armstrong & German 2000, Lane, 
Creech-Eakman & Nordgren 2002). 

The latter method is more direct and less model dependent but 
up to now it has been applied only to a limited number of stars. 
On the other hand, the different versions of the BW technique 
can be applied to relatively large Cepheid samples but require 
both accurate photometric and radial velocity data. 
A powerful modification of the BW technique is the so called 
CORS method (Caccin et al. 1981), which has the advantage of 
taking into account the whole light curve rather than selecting 
phase points at the same color (as in the classical BW imple- 
mentation), but relies on the adoption of an accurate Surface 
Brightness (SB) calibration. 

Originally Sollazzo et al., (1981) adopted the empirical SB 
photometric calibration in the Walraven system provided by 
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Pel (1978). More recently Ripepi et al. (1997) modified the 
method, by adopting the empirical calibration of the reduced 
surface brightness Fy as a function of (V - R) provided by 
Barnes, Evans & Parson (1976). This modified version of 
the CORS method was tested, for different colors selections, 
through the application to synthetic light and radial velocity 
curves based on nonlinear convective pulsation models (Ripepi 
et al. 2000). 

The recent release of new Cepheid data in the Stromgren pho- 
tometric system (Arellano-Ferro et al. 1998), and the known 
sensitivity of intermediate band colors to stellar physical pa- 
rameters (e.g. gravity and effective temperature) suggested us 
to investigate the possibility of extending the CORS method to 
the Stromgren filters. 

To this purpose we have derived, in this system, a SB cali- 
bration based on model atmosphere tabulations. In this paper 
we present a modified version of the CORS method based on 
this new cahbration and the application to a sample of Galactic 
Cepheids. 

The organization of the paper is the following: in Sect. 2 we 
summarize the assumptions and the philosophy of the tradi- 
tional CORS method; in Sect. 3 we introduce the modified 
CORS method based on the new SB calibration; in Sect. 4 we 
test the new method by means of pulsation models; in Sect. 
5 we apply the method to a sample of Galactic Cepheids and 
present the comparison with the literature; in Sect. 6 our final 
results concerning the PR relation for Galactic Cepheids are 
shown and the theoretical fit of observed light and radial veloc- 
ity variations for Cepheid Y Oph is used as an additional check. 
Some final remarks close the paper. 

2. The original CORS method 

In this section we briefly outhne the assumptions and the main 
features of the CORS method in order to better understand what 
follows. 

The CORS mefliod (Caccin et al., 1981) starts from the defini- 
tion of the surface brightness: 

Sv + 5-\oga (1) 

where a is the angular diameter of the star. For a variable star, 
Eq.Q is valid for the whole pulsational cycle, so that differen- 
tiating it with respect to the phase (0), multiplying by a color 
index (C,j) and integrating over the whole cycle, one obtains: 

q ^ Inl^Ro-kP J v(0') d(f>'^ C-jdcf) 

-B + AB = (2) 
where q ^ 

B= r Cij(cf>)m'y(<p)d4> (3) 
Jo 

AB= r CtjWS'yWdcf> (4) 
Jo 

where P is the period, v the radial velocity and k is the radial 
velocity projection factor which relates radial to pulsation ve- 
locity (R'((/>) - -k-P- v(0)). The typical value for ^ is 1 .36 (see 



discussion in Ripepi et al. 1997). 

From Eq.(|3 we can evaluate the radius at an arbitrary phase 
00 (usually taken at the minimum of the radial velocity curve), 
whereas the mean radius is obtained by integrating twice the 
radial velocity curve. 

The B term can be easily calculated from observed light and 
color curves. On the contrary, the AB term, which includes the 
Surface Brightness, is not directly observable. By neglecting it 
in Eq. Q, we obtain the pure Baade-Wesselink method (see 
Caccin et al. 1981). However, Sollazzo et al. (1981) and Ripepi 
et al. (1997, hereinafter RBMR) demonstrated that the inclu- 
sion of AB improves the accuracy of radius estimates, provided 
that 5 V is evaluated at each pulsation phase. 

3. An Improvement of the CORS method based on 
the Stromgren photometry 

3.1. Evaluation of the AB term 

As outlined in the previous section, the inclusion of AB im- 
proves the accuracy of radius estimates. In this section we 
present a new good approximation for this term. 

As discussed by Onnembo et. (1985), if the quasi-static ap- 
proximation (QSA) ' is assumed for Cepheid atmospheres, any 
photometric quantity can be expressed as a function of effective 
temperature and gravity (Te/f, geff)\ Then we can write: 

Sv ^ Sv(Teff,geff) (5) 

and 

Cij^Cij(Teff,geff) (6) 
Ckl = Cu(Teff,geff) 

where 5v is the surface brightness in the visual band, and 
Cij = (;«; - nij), Ckl - (nik - mi) are two arbitrary colors. 
If the last two equations are invertible, i.e the Jacobian 
J{Teff,geff\Cij,Chd + 0, then we can invert Eqs. (|6j, obtain- 
ing: 

T^'^fi = T^eiMij, Ckl) (7) 
8eff = 8eff(Cij, Ckl) 

and hence 

Sv ^ Sv(cij,Cki) (8) 

We point out that Eqs. (|6j do not admit a general solution over 
the whole parameter space, given that the same colors can be 
obtained for different pairs of Tgf f and gc//. This notwithstand- 
ing it is possible to find a local solution, that is a solution valid 
only in the parameter space defined by Cepheids. 

A potentially advantageous choice for the colors Cjj, Cki 
could be represented by the Stromgren reddening free indexes 

' We recall that a sufficient condition for the validity of the QSA is 
that the atmosphere of the pulsating star can be described, at any time, 
by a classical hydrostatic, plane parallel model in radiative/convective 
equilibrium and local-thermodynamic-equilibrium (LTE), identified 
by the effective temperature T^// and by the effective gravity g^ff = 

GM , £r 
«2 A2 
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Fig. 1. A grid of lines at constant Tejf and constant log g in the 
theoretical [oti], [ci] plane. 

[mi] and [ci] defined as follows (Crawford & Mandwewala 
1976): 

[mi]^mi+0.33{b-y) (9) 

[ci] =ci- 0.16(^7 -y) (10) 

where the coefficients 0.33 and -0.16 are suitable for F 
supergiant stars (Gray 1991). 

In this context, an interesting possibility to derive Eqs. 
0, consists in using grids of theoretical colors, calculated by 
means of model atmospheres, to obtain T^/f and g^jj as a func- 
tion of [oti] and [cj]. 

To verify this possibility, we have adopted the grids of theoret- 
ical colors by Castelli, Gratton & Kurucz (1997a,b, hereinafter 
C97a,b). The theoretical grid (constant T^/f and \ogg) in the 
[mi], [ci] plane is shown in Fig.^ This figure suggests that for 
0.50 < logg < 3.50 and 5000/: < r,// < 7000/: there is a one- 
to-one correspondence between a point in the [oti], [ci] plane 
and the corresponding T^f / and log g values. Thus, in principle, 
in this color range it is possible to invert Eqs. (|6j and, in turn, 
to derive an expression for Eq.(|^. 

Before proceeding, it is important to verify if the location 
of theoretical grids in the [mi], [ci] plane is consistent with the 
one occupied by real Cepheid data. To this aim, we have over- 
plotted on Fig.^the color-color [mi], [ci] loop for all the stars 
in our sample (see section 5.1). Figure |2 shows the resulting 
comparison for three stars (FF Aql, FN Aql, U Aql) of our 
sample, characterized by short and intermediate periods. The 
figure shows that the color-color loops for the selected stars 
are completely included in the theoretical grids. The same test 
has been performed for the three stars W Sgr, WZ Sgr, S V Vul, 
which have longer periods. As showed by Fig.|3] in this case the 
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Fig. 2. Theoretical grids in the [mi], [ci] plane, compared with 
the empirical loops of three stars with period shorter than 10 d. 
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Fig. 3. The same as in Fig. 2 but for three Cepheids with period 
longer than 15 d. 



color-color loops for the selected Cepheids lies outside (at low 
gravity and low effective temperature) of the region covered by 
theoretical grids. This unexpected result has to be taken into 
account when applying our method to such long period stars. 
An explanation for this model limitation is beyond the scopes 
of present paper but it is an important issues worth to be ad- 
dressed in a future work. 
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Fig. 4. log Teff as a function of [mi],[ci] obtained by a 4 
degree polynomial fit of theoretical data (crosses) by C97a,b. 

Fig. 5. The same as in Fig. 4 but for log g. 
3.1 .1 . Derivation of 5\/ 

As shown in the previous section, the first step in the con- 
struction of a new, based on Stromgren photometry, version 
of the CORS method, is the formulation of 5 v in the form 
of Eq.ljSji. To this aim, we have first to find relations in the 

form of Eqs.(0 starting from theoretical grids. This has been 

th 

achieved by means of a 4 degree polynomial fit to the 
theoretical grid (least square fit). The results of the polynomial 
approximations are reported in Appendix A, and showed in 
Figs-Sand|5lfor log T^jf and log g^jf respectively. The rms of 
the fits are:0.0018 dex for log T^ff and 0.1 for logg^ff. 

On the basis of the calculated relations we are now in the 
position to estimate the surface brightness from the expression: 

Sv = const. - IQlogT, ff - BC(Teff,geff) (11) 

where BC is the Bolometric CoiTection, derived as a function 
of effective temperature and gravity {BC - BC{Tgff, ggff)). 
This function can be easily obtained using again a 4^^ degree 
polynomial fit to the theoretical grids. The resulting equation is 
also reported in Appendix A and shown in Fig|6l Note that the 
rms is only 0.003 mag. 

The procedure outlined above allowed us to achieve 
our goal, i.e. to derive an analytic, although approximated, 
expression of Sy = 5'i/([mi], [ci]), and, in turn, to calculate 
the AB term, which allows us to apply the CORS method, 
in its more general formulation, for data in the Stromgren 
photometric system. 



Fig. 6. BC as a function of logT;,, \ogg obtained by a 4 degree 
polynomial fit of theoretical data (crosses) by C97a,b. 



2. pulsation models provide a set of equally well sampled 
curves covering wide period and effective temperature 
ranges; 

3. for pulsation models the radius and intrinsic luminosity in- 
formation are available. 

In the following we describe in detail how the models have 
been used to test the method. 

4.1. The synthetic curves 

To test both the accuracy and the consistency of the new CORS 
method we have adopted the pulsation observables predicted 
by hydrodynamical models of classical Cepheids. A detailed 
discussion on the physical assumptions adopted to calculate 
these models can be found in Bono, Marconi & Stellingwerf 
(1999, BMS99) and Bono, Castellani, Marconi (2002). Among 
the different sequences of nonlinear models computed by 
BMS99 we have selected canonical models^ at solar chemical 
composition {Y = 0.28, Y = 0.02) and stellar masses ranging 
from 5 to 9 Mq. At fixed stellar mass, we generally chose 
three models which are located in the middle of the instability 
strip as well as close to the blue and the red edge. The period 
of selected models roughly ranges from 3.5 to 62 days. The 
input parameters (M,Teff, L), the computed radius and the 
pulsational period are summarized in Tabled 



Table 1. Physical properties of the selected Cepheid models 



model 


Mass 


Luminosity 


Teff 


Radius 


Period 




Mo 


logL/Lo 


K 


Re 


Days 


modi 


5 


3.07 


5800 


34.3 


3.5091 


mod2 


5 


3.07 


5600 


36.7 


3.9412 


mod3 


6.25 


3.42 


5400 


58.6 


7.5842 


mod4 


6.25 


3.42 


5100 


65.2 


8.7060 


mod5 


7 


3.65 


5300 


81.0 


12.0904 


mod6 


7 


3.65 


5000 


91.3 


14.7582 


mod7 


7 


3.65 


4800 


97.3 


16.8322 


mods 


9 


4.0 


4900 


185.7 


46.5687 


mod9 


9 


4.0 


4700 


204.1 


53.9378 


mod 10 


9 


4.0 


4500 


220.6 


62.3650 



4. Test of the revised CORS method 

To verify the accuracy of the new approximation, we have 
applied the CORS method (with the new Sy calibration) to 
synthetic light, color and radial velocity curves predicted by 
Cepheid full amplitude, nonlinear, convective models. The ad- 
vantages of testing the method by means of pulsation model 
curves, rather than empirical data, are the following: 

1. for pulsation models, bolometric light curves are trans- 
formed into the observational bands using the same model 
atmosphere tabulations (C97a,b) we have used to derive the 
surface brightness calibration; 



Theoretical observables have been transformed into the ob- 
servational plane by adopting the bolometric corrections (BC) 
and the color-temperature relations by C97a,b. We assumed 
Mboi(Q) = 4.62 mag and adopted atmosphere models com- 
puted by neglecting the core overshooting and for a fixed value 
of the microturbolence velocity ^ - 2kms^^ . 
For each pulsation model we have derived the V light curve, 
three color curves, namely in (b-y), [mi] and [ci] and the radial 
velocity curve. Figure shows these curves for models with 
M = 9M0, Teff = A5QQK, and M = IM^, T^ff = 53QQK. 

^ Canonical models are the ones constructed by adopting a mass- 
luminosity relation based on evolutionary computations which ne- 
glect convective core overshooting during hydrogen burning phases 
(Castellani, ChiefR & Straniero 1992). 
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Fig. 7. Variations along a full pulsation cycle of magnitude, 
color and radial velocity for models with M = 9Mq and Te/f = 
4500/: (left) and M = 7Mq and T^ff = 5300/: (right). 

4.2. Application to theoretical models 

As a first test, we have applied the method to "perfect" light, 
color and radial velocity curves, i.e. the curves directly ob- 
tained from the models which, of course, do not show random 
errors. This test allows us to verify whether or not our calibra- 
tion of the surface brightness is intrinsically correct. If this is 
the case we do not expect a large discrepancy between the cal- 
culated radii and the "true" (theoretical) ones. 
To perform this comparison we have evaluated the radius of 
each model in the two different approximations: 

- Without the AB term (Baade-Wesselink approximation); 

- With the AB estimated as described in section 3. 

The results of this test are summarized in the first four 
columns of Table |2 We note that, the application to models 
with long period is justified by the fact that the bolometric light 
curves of pulsation models have been transformed into the ob- 
servational bands by using the same model atmospheres by 
C97a,b. Therefore, at variance with the observed long period 
Cepheids, the model color-color loops are always consistent 
with C97a,b's grids. 

The same result is also shown in Fig. |S| which suggests 
that our calibration of the surface brightness is intrinsically 
coiTect with the discrepancy between "computed" and "theo- 
retical" Cepheid radii being around 1 %. We also note that the 
inclusion of the AB term only slightly improves the agreement 
with predicted radii, producing a small reduction of the scatter 
around the mean (see labelled values in Fig.|S}. 

4.3. Stability tests 

Once verified that the new calibration of the surface brightness 
is intrinsically correct, we are in the position to estimate the sta- 
tistical error in the radius determination due to uncertainties in 
the measurements, i.e. to test the sensitivity of our technique to 
the accuracy of observed light, color and radial velocity curves. 
To this aim, we have transformed the synthetic curves for mod- 
els in Table [Oby adopting the following steps: 



Table 2. Radii estimated with the two different approximations 
, i.e with and without the AB term (respectively column 2,3), 
using model predictions (column 4), and from modified syn- 
thetic curves in the case of good (column 5) and fair (column 
6) data respectively. 



Models R„iii,oi«AB 


Ras 


Rrfo 




Ras 


R 


R 


R 


R 


R 




Ro 


Ro 


Ro 


Ro 



modi 
mod2 
mod3 
mod4 
mod5 
mod6 
mod7 
mods 
mod9 
modlO 



34.7 

36.7 

59.1 

64.6 

82.1 

92.4 

98.6 

190.3 

206.9 

223.1 



34.8 
36.6 
59.0 
65.7 
82.2 
92.3 
98.5 



34.3 
36.7 
58.6 
65.2 
81.0 
91.3 
97.3 



37.5 
36.7 
59.5 
63.1 
81.3 
90.6 
99.4 



36.3 
38.6 
60.9 
73.1 
78.0 
94.0 
99.6 



188.8 185.7 185.8 1534 
206.5 204.1 202.0 203.1 
223.4 220.6 187.7 177.0 



m=1.010 (T=0.009 



m=1.010 (7=0.006 



Fig. 8. Ratio between "computed" and "theoretical" radii as 
a function of the logarithmic period. The bottom panel dis- 
plays the radius evaluation based on the revised CORS method, 
while the top one the radius evaluations based on the pure BW 
method. 



1 . we have reduced the number of phase points to the typical 
value for observations, i.e. ~ 30 - 35; 

2. we have extracted phases randomly; 

3. we have added Gaussian errors, considering the cases: 

a) good data (cry = O.Olmag, cr[,„^-\ — 0.04mag, cr[^.,] = 
Q.OAmag, ctrv - 0.25A:to"', number of phase points=35); 

b) fair data {crv - 0.04-mag, cr[„,j] - O.OSmag, o-[c,] - 
O.OSmag, ctrv - 0.5A:m"', number of points=30). 

As an example, Fig.|9]shows the synthetic light, color and 
radial velocity curves obtained for model mod4 (see TableQ. 
Left panels show the pure model curves, whereas middle and 
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Fig. 11. Ratio between the "computed" and "theoretical" radii 
as a function of the shift A<p for the model (mod2) with M = 5 
Mo, Teff = 5600/:. 



Fig. 9. Light, color and radial velocity curves for the model 
(mod4) with M = 6.25Mq and T,ff = 5100/:. The left pan- 
els show the pure model curves; the middle and right ones the 
modified curves in the case of good and fair data respectively. 



m=1.02 (T = 0.05 



I I I I I I I I I 
L m=0.99 (T = 0.06 



m=1.02 (7 = 0.10 



. L m = 0.99 (T = 0.09 



1 1.5 

logP 



logP 



Fig. 10. Ratio between "computed" and "theoretical" radii as 
a function of the logarithmic period. In the left panels (a) the 
radius is computed using modified curves that simulate good 
quality, adopting either the revised CORS method (bottom) or 
the pure BW one (top). In the right panels (b) similar plots are 
shown for radii computed from modified curves that simulate 
fair quality data. 



right panels show the modified curves in the two cases a) and 
b) respectively. 

We then applied the new CORS method to the whole set of 
modified synthetic curves. The resulting radii are reported in 
the last two columns of Table|2 whereas in Fig.fTOjwe plot the 



ratio between "computed" and "theoretical" radii as a function 
of the logarithmic period. In particular, top and bottom pan- 
els show the CORS solutions without and with AB respectively 
for the cases a) (left) and b) (right). Figure Q^seems to show 
that including the AB term does not improve much the results, 
on the contrary, the scatter in case a) worsen. However, this 
disagreeable occurrence is only apparent. In fact, by exclud- 
ing "modS" and "mod 10" that show a very peculiar morphol- 
ogy (a sharp bump) of the [ci] color curve (see figure Fig0for 
"mod 10"), which makes difficult the fit and the AB calculation, 
the average uncertainties on the radius estimation fall to 3.5% 
and 5% in cases a) and b) respectively (including the AB term). 
These numbers represent a useful lower limit for the eiTor asso- 
ciated with radius determination obtained by the CORS version 
developed in this paper. 

Another possible source of uncertainty in the radius deter- 
mination is the misalignment between hght or color curves, and 
the radial velocity one. In fact, the photometric and radial ve- 
locity data are rarely collected simultaneously. This occuiTence 
could introduce a shift in phase A<p between the two different 
data set. To verify the importance of this shift on the radius de- 
termination we have introduced an artificial phase shift (start- 
ing steps of A(f) = 0.01 up to A0 = 0.1) in the synthetic radial 
velocity curves with respect to the photometric ones. 

The result of such a test is shown in Fig.^] where the ra- 
tio between "computed" and "theoretical" radii is plotted as a 
function of the phase shift A(p for the model mod2. The "com- 
puted" radius is larger than the "theoretical" one for < Acp < 
0.1, whereas it becomes smaller beyond 0.1 with the effect in- 
creasing with the phase shift. Eventually, our program does not 
converge for A0 > 0.3. However, such large phase shifts are 
never reached when we deal with actual data. In fact, when the 
temporal distance between photometric and radial velocity data 
is rather long (~ 1000 cycles ), for a typical Cepheid with pe- 
riod P ~ 10 + 5 ■ 10"^ d, the resulting phase shift is A0 ~ 0.05. 
As shown in Fig. this value of A<p is sufficient to generate 
a systematic error on the derived radius of about 17%. This oc- 
currence confirms quantitatively the need to use photometric 
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and radial velocity data as close as possible in time, or other- 
wise to correct this shift. 

5. Application of the method 

Having tested the capabilities of the new version of the CORS 
method (see previous section), we are ready to apply it to actual 
data. In the following we discuss the application to a sample of 
52 Galactic Cepheids. 

5.1. The sample 

We searched the literature for an homogeneous sample of 
Galactic Cepheids with photometric data in the Stromgren sys- 
tem and we selected the 31 pulsators analysed by Arellano 
Ferro et al. (1998) (AFGR hereafter) supplemented by data for 
other 21 objects from the papers by Feltz & McNamara (1980) 
(FM hereafter) and Eggen (1985, 1996) (Eg hereafter). 

Concerning radial velocity data, several large and homo- 
geneous datasets are available in the literature. In particular, 
we have used the catalogues from: Evans (1980) (E hereafter), 
Gieren (1981) (G hereafter), Barnes et al. (1987, 1988), (BMS 
hereafter), Coulson et al. (1985) (CCG hereafter), Wilson et 
al. (1989) (WCB hereafter), Metzger et al. (1993) (MCMS 
hereafter), Bersier et al. (1994) (BBMD hereafter), Gorynya et 
al. (1998) (GSSRGA hereafter). Kiss (1998) (K hereafter) and 
Imbert (1999) (I hereafter). 

In general we considered only Cepheids with good pho- 
tometric and radial velocity curves, i.e. more than 15 phase 
points and reasonable precision. When measurements by dif- 
ferent authors are available for the same stars we have chosen 
data which have: 1) the largest phase points number and the 
best precision; 2) the lowest temporal separation with respect 
to the photometric data. When we have distant data sets, that 
is differences between radial velocity and photometric curves 
larger than ~ 1000 cycles, (6 stars, see Table O we correct 
the misalignment between this two curves by using different 
epochs. 

In few cases we have merged data from different authors, in 
order to obtain more sampled radial velocity curves (see Table 
|3}- Our final selected dataset is summarized in Table|3l 

5.2. Peculiar stars 

Some Cepheids in our sample need to be discussed individu- 
ally: 

- For the stars TT Aql, X Cyg, T Mon, RU Set, Y Set, WZ Sgr 
and SV Vul, we can calculate the radius only in the classi- 
cal Baade-Wesselink approximation (without the AB term), 
because for these stars the loop in the [mi], [ci] plane does 
not lie completely within the theoretical grids by C97a,b 
(see paragraph 3.1); 

- We have excluded from our sample the Cepheids GI Car, 
SU Cas, Y Lac, SZ Tau, X Vul for the poor data quality. 

- Stars DT Cyg, FF Aql, BG Cru and V440 Per are first over- 
tone pulsators (Antonello, Poretti & Reduzzi, 1990). 



- Stars TU Cas, V367 Set and BQ Ser are double mode 
Cepheids (Pardo & Poretti, 1997 and references therein). 

- Stars U Aql, FF Aql, V 496 Aql, 77 Aql, RX Aur, SU Cas, 
6 Cep, GI Car, SU Cyg, VZ Cyg, ^ Gem, X Lac, Y Lac, Z 
Lac, BG Lac, T Mon, Y Oph, BE Oph, AW Per, Y Set, 
S Sge, W Sgr, Y Sgr, WZ Sgr, AP Sgr, V 350 Sgr, SZ 
Tau, T Vul, U Vul and SV Vul are members of binary sys- 
tems. Note that this list is the result of a detailed analy- 
sis of Szabados's selection. (Szabados, 2003 and references 
therein). 

5.3. Comparison with the literature 

In table|4]we report the radius obtained for all the Cepheids in 
our sample, and, for comparison purposes, we report the liter- 
ature results renormalized to our projection factor k. In partic- 
ular, from left to right we report for each Cepheid: the name; 
the period; the CORS radii obtained without and with the AB 
term (RmthoutAB and Rab); the radii obtained by Arellano Ferro 
& Rosenzweing (2000) Rafr, Laney & Stobie (1995) Ris or 
Gieren et al. (1998) Rcfc, Ripepi et al. (1997) Rrbmr. and 
Moffett & Barnes (1987) Rbm- 

In Fig. [21 we have compared our results with those of authors 
reported in the table 0] As a result we find that our radii are 
on average larger than the radii obtained by Arellano FeiTO & 
Rosenzweing (2000) and Moffett & Barnes (1987), while they 
are slightly smaller than Ripepi et al. (1997) ones (see labels 
in Fig. I12> . The scatter in this comparisons is rather large (~ 
30%); this occurrence could be due to: 1) the inclusion of bi- 
nary stars in the comparison (different methods uses different 
colours and then binarity could affect differentially the various 
determinations); 2) the use of optical colors. Verifying hypoth- 
esis 1) by comparing only the non-binary stars does not make 
sense because of the small number of "bona fide" single stars 
in our sample. 

Concerning point 2), as well known in the literature (Laney 
& Stobie 1995, Gieren, Foque & Gomez 1997), radii obtained 
from NIR data are more precise than the ones obtained by using 
optical colors. Unfortunately, most works with radii determina- 
tions from NIR data deal with southern Cepheids and we have 
very few stars in common (see column 7 of Table 4 for the 
comparison) 

Recently it has become possible to derive accurate stellar 
radii from the angular diameters measured with interferometric 
techniques (see Nordgren et al. 2000 and Lane et al. 2002) com- 
bined with Hypparcos parallaxes. The comparison with these 
measurements represents a useful test for our results. In Table 
|5]we compare with our results the radii obtained by Nordgren 
al. (2000) and Lane et al. (2002), with the quoted interfero- 
metric techniques, for the stars in common with our sample, 
namely Aql, 6 Cep and ( Gem. In particular, from left to 
right, we report: the name of the star; the radius {Rj) obtained 
by means of interferometric techniques {6 Cep from Nordgren 
et al. 2000 and //Aql and ( Gem from Lane et al. 2002); our 
radius in the two different approximations, without {RwnhouthB) 
and with (RKirhAB) the AB term respectively. The interferomet- 
ric radii, reported in Table |5] are corrected for the different k 
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Table 3. The selected Cepheid sample 



Cepheid 


Period 


Photometry 


Kad. Vet. 


i\(rn — KV) 


Cepheid 


Period 


Photometry 


Kacl. Vet. 


isyrn — KV) 




days 


source 


source 


cycles 




days 


source 


source 


cycles 


FF Aql 


4.470840 


AFGR 


BMS 


1142 


YOph 


17.126780 


AFGR 


GSSRGA 


46 


FM Aql 


6.114334 


AFGR 


GSSRGA 


89 


AW Per 


6.463589 


FM 


GSSRGA 


422 


FN Aql 


9.481614 


AFGR 


GSSRGA 


56 


V440 Per 


7.572712 


Eg 


GSSRGA 


131 


77 Aql 


7.176779 


FM 


BMS 


413 


CM Set 


3.916977 


AFGR 


M CMS, B 


463 


TT Aql 


13.75551 


AFGR 


GSSRGA 


81 


EVSct 


3.090998 


AFGR 


BBMD, MCMS 


586 


UAql 


7.024100 


AFGR 


BMS 


674 


RUSct 


19.70062 


AFGR 


MCMS 


92 


V496 Aql 


6.807164 


AFGR 


G 


643 


SS Set 


3.671280 


AFGR 


BMS, G 


1194 


V600 Aql 


7.238748 


AFGR 


GSSRGA 


156 


V367 Set 


6.29307 


AFGR 


MCMS 


317 


RT Am- 


3.728561 


FM 


GSSRGA 


216 


YSct 


10.341650 


AFGR 


BMS 


421 


RX Aur 


11.623537 


FM 


1 


271 


BQ Ser 


4.316700 


AFGR 


GSSRGA 


926 


GlCar 


4.431035 


Eg 


G 


- 


S Sge 


8.382044 


FM 


GSSRGA 


209 


SUCas 


1.949322 


FM 


GSSRGA 


- 


AP Sgr 


5.0574269 


AFGR 


BMS, G 


864 


TU Cas 


2.139298 


FM 


GSSRGA 


- 


BB Sgr 


6.637115 


AFGR 


GSSRGA 


166 


<5Cep 


5.366316 


FM 


BBMD 


393 


USgr 


6.745363 


AFGR 


GSSRGA 


81 


DT Cyg 


2.499086 


AFGR 


GSSRGA 


2203 


V350 Sgr 


5.154557 


AFGR 


GSSRGA 


157 


SUCyg 


3.845733 


AFGR 


BMS 


1432 


XSgr 


7.012630 


AFGR 


WCB 


759 


XCyg 


16.386130 


AFGR 


BBMD 


244 


YSgr 


5.773400 


AFGR 


BMS, WCB 


885 


VZCyg 


4.864504 


FM 


BBMD 


50 


YZSgr 


9.553606 


AFGR 


BMS 


419 


BGCru 


3.342503 


Eg 


E 


943 


WSgr 


7.595080 


AFGR 


BBMD 


542 


WGem 


7.913960 


FM 


I 


257 


WZSgr 


21.849708 


AFGR 


GSSRGA 


50 


f Gem 


10.150780 


FM 


BBMD 


179 


SZ Tau 


3.148727 


FM 


GSSRGA&- 




BGLac 


5.331938 


FM 


I 


114 


SV Vul 


44.999660 


AFGR 


1 


103 


YLac 


4.323776 


FM 


I 


46 


T Vul 


4.435532 


FM 


BBMD 


530 


Z Lac 


10.88554 


FM 


BMS,GSSRGA 


273 


U Vul 


7.990736 


AFGR 


GSSRGA 


142 


XLac 


5.444990 


FM 


BBMD 


1.9 


X Vul 


6.319562 


FM 


GSSRGA 




TMon 


27.024649 


FM 


GSSRGA 














BFOph 


4.067695 


AFGR 


BMS, G 


1077 













Photometry sources: Arellano Ferro et al. (AFGR,1998); Feltz & McNamara (FM,1980) and Eggen (Eg,1985,1996). 

Radial velocity sources: Evans (E,1980); Gieren (G,1981); Barnes et al. (BMS,1987, 1988); Coulson et al. (CCG,1985); Wilson et al. 
(WCB, 1989); Metzger et al. (MCMS,1993); Bersier et al. (BBMD,1994); Gorynya et al. (GSSRGA, 1998); Kiss (K,1998) and Imbert (1,1999). 
A(P/i - RV): number of pulsational cycles between photometric and radial velocity curves. 



Table 5. Comparison with radii measured with interferometric tech- 
niques (see text for details). 



Star 


Ri 


^YvillwutAB 


RwithAB 


77 Aql 


58.8 ±7.6 


48.2 


52.6 


(5Cep 


471 


41.0 


47.2 


(TGern 


63.4 ± 7.2 


64.4 


75.0 



projection factors. So we obtain the radii 61.8 x 1.36/1.43= 
58.8 /?o, 45x1. 36/1. 31=47/^0 and 66.7 x 1.36/1.43=63.4 /?o for 
Aql, 6 Cep and ( Gem, respectively. We notice that the agree- 
ment between our results and the interferometric ones is very 
good for 77 Aql and 6 Cep, whereas some discrepancy is found 
for ^ Gem . However, by assuming a typical error of ~ 10% (see 
discussion in Sect.7) on the radii found, we conclude that our 
results are globally consistent with the interferometric ones. 



6. The Period-Radius relation 

After having tested the method and the consistency of our re- 
sults with previous determinations in the literature, we are in 
the position to derive a PR relation based on the new derived 
radii. 

The PR relations for all the stars in our sample for which the 
program reached the convergence are shown in Fig.^J where 
the the top and bottom panels illustrate the case without and 
with AB respectively. A least square fit to the data (solid lines 
in Fig. I13> leads to the following Period-Radius relations in 
the case without the AB term (Eq. ll2t and with the AB term 
(Eq.in}: 

logT? = (1.18 + 0.05) + (0.67 + 0.06) log P r.m.i = 0.09 (12) 

log /? = ( 1.22 + 0.06) + (0.67 + 0.08) log P r.m.i = 0.08 (13) 

An inspection of the figure suggests that the presence of 
first overtone pulsators (filled triangles) could affect the derived 
PR relations because, at fixed period, they are brighter and, in 
turn, have larger radii than fundamental pulsators. 
Similarly, also Cepheids belonging to binary systems could 
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Cepheid 


Binary 


Period 
days 


R_ 

Rq 


Ras 

_R_ 

Rq 


^AFR 
_R_ 

Rq 


^LS 1 ^GFG 
_R_ 

Rq 


^RBMR. 
_R_ 

Rq 


_R_ 

Rq 


Anl* 
rr rtCJi 


\j 




AQ 8 


'>A 1 
J'4. 1 


'>9 n 

JZ.U 














SA 




'>Q ^ 




AH A 
DU.O 


SI Al 
J 1 .01 


rIN Aqi 








8A ^ 










77 Aql 


R 




48 9 


'^9 A 
JZ.O 


A7 ^ 
/ . J 




'^A 7 
JO. / 


S9 7A 
JZ. /O 


TT An^ 
1 1 /\C[i 




ID. ID JO i 


04. :> 




7A 1 
/0. 1 




Q^ ^ 
VJ. J 


QA ^ 1 
y4. J 1 


TT Anl 


kJ 


7 nozi 1 nn 

/ .UZ'4-iUU 


A^ A 


A7 n 


A8 Q 




Al ^ 
01. J 


S9 AS 
JZ.4J 


VdQA Anl 


R 


A 807 1 fiA 


Al 7 


Al 9 
D 1 .Z 


^8 1 
JO. 1 




'lA A 
JO. 4 


Al AQ 
4J.07 


\7f.r\r\ Anl 
V OUU /\C[1 






J J. 4 


8 n 1 

U.l 


'^A ^ 
JO. J 








PV Aiir 
r\/\ /\Ui 


R- 


1 1 A9'^'^^7 


S'^ 7 


81 n 
1 .U 








Al 17 
01.1/ 


RT Aiir 




J. 1 Zo JDi 


97 Q 

L 1 .y 


■^n A 

JV.O 


^A 9 
jO.Z 




A9 9 
4Z.Z 


^S 7 A 
J J. /4 




V 


^ QAAa 1 A 
J. jDOj 10 


A\ n 


Al 9 
4 / .Z 


AA 9 
40. Z 




^9 8 
jZ.o 


^Q aA 






Z.H-77UoD 


AA A 


^8 7 


■^9 ^ 
JZ. J 




A9 Q 
4Z.7 






n 






AO 8 


A8 Q 
40. 7 




Al ^ 
01. J 


S9 AS 
jZ. 4J 






1 A ^SA 1 


1 OA 1 
1U4. 1 




AQ ^ 
O7. J 




QA 9 
VO.Z 


1 1 A 9Q 

1 14. zy 




r\ 
yj 


A 9.f\A'>C\A 


A'^ '\ 


J J. 1 


8 




A7 
4 / . J 


^S S7 

J J. J / 






J. jH-ZjUj 


98 Q 
Lo.y 


■^9 7 
^Z. / 












V 


1 n 1 '>n7sn 
lu. ijy) 1 oU 


AA A 
04. 4 


7^ n 

/ J.U 


J J. 1 




8A 9 
oO.Z 


A9 SS 

OZ. J J 


W Ocill 




7 Ql ^QAO 
/ .71 jyou 


QA Q 


8A A 
oO.'4 


A8 S 
4o. J 




AH 7 
OU. / 


A8 AQ 




u 
u 


J. J J lyjo 


Al n 


A8 Q 






8 

JJ.O 


A'^ AQ 
4 J. 47 


V T n<- 


D 




Al 'X 


AA Q 


■^A 7 
jO. / 




71 9 
/ 1 .Z 


Al 07 


7 T dr- 




lU.oo J JH- 


A7 9 
D / .Z 


AQ A 


8^ 1 
J. 1 




8A 8 
oO.o 


AA AO 
00. 4U 


Dr vjpn 


D 




Dy.y 


AA 8 




■^A 1 1 / a'' 

JO. 11/ J J.o 


'?9 A 
jZ.4 


78 
J J. /o 


1 L/pn 


D 


1 7 1 9A7Cn 
1 / . IZO / oU 


QA ^ 

yo. J 


Q^ A 




yj.LL 


119 9 
1 IZ.Z 




AW Pfir 
rWV r CI 






9Q 7 
ZV. / 


^7 S 


■^A 7 
j4. / 




'\A 8 
J4.0 


AS AH 
4J.4U 


V'+^U rcl 


R ■ 


7 ^797 1 9 
/ . J /Z / IZ 


A9 Q 

oz.y 


7A a 

/O.O 














Q Q 1 AQ77 


^8 A 


A 

JO.O 


9Q A 
Z7.O 








S-j V ijCL 


R- 






A9 n 




Al A9«/ "^9 
4/.OZ / jZ.j 


■^A A 
jO.O 








^ A71 980 
j.O / IZoU 


A 


AA 9 


in A 






97 71 
Z / . / 1 


PTT ^nt 




1 Q 7nnA9 


1 09 Q 






1 90 A'>^ 
1ZU.4J 






V *i/-t 


K 

D 


iU. JH- 1 DJU 


AO Q 

ou.v 








71 '^ 
/ 1 . J 


80 AS 
oU.4j 


Ap c„- 


R ■ 


c n^7d9AQ 


An A 


A9 ^ 
OZ. J 


^A A 
jO.O 




jU.U 




5 oge 


\j 


0. JOZU'4'+ 


7^^ 1 

/ J. 1 


7Q 1 


AA A 
04. 








RR tlrrr 




A A'^71 1 
0.0^ / 1 1 J 


Al A 


Al ^ 
01 .J 




■^7 /A A A^' 
J 1 .Jy /44.4 


J J. J 


AO A9 
4U.4Z 


u ^gr 




0. / 4J jOj 


7A A 
/O.O 


89 ^ 
oZ. J 




'^9 ^n^ / A8 a'' 
JZ.jU / 4o.o 




S7 Al 
J / .01 


V jjju ogr 


r\ 
yJ 


S 1 '^A'>S7 
J. 1 JH-J J / 


^9 1 
jZ. 1 


^9 7 


AO 9 
4U.Z 




AA n 

40. U 


AQ SS 
47. JJ 


A 5gr 




7 ni 9A^n 

/ .UIZDjU 


S9 


'\Q 8 


A8 S 
Oo. J 




'>8 A 
Jo.O 


A7 7A 
4 / . / 


Y Sgr 


D 


Q TV Q /inn 


A 

4/A 


JO.Z 






01.4 




YZ Sgr 




9.553606 


93.6 


95.8 






136.0 




WSgr 





7.595080 


58.6 


62.6 


62.5 




58.6 


60.77 


WZ Sgr 


B 


21.849708 


120.5 






121.43^^ /122.2'-' 




151.20 


SV Vul 


B 


44.999660 


213.0 




135.6 


241.627 250.7^ 




189.5 


T Vul 


b 


4.435532 


31.8 


37.1 


42.3 




48.8 


36.54 


U Vul 





7.990736 


75.0 


84.4 


57.9 




61.1 


54.34 



An asterisk at top-right of Cepheid name means that it pulsates as a First Overtone. 

In the second column: B - spectroscopic binary, B: - spectroscopic binary but confirmation needed; b - Photometric companion, physical relation 
should be investigated; O - spectroscopic binary with known orbit; V visual binary (see Szabados, 2003 for details). 

Radius sources: Rafr Arellano Ferro & Rosenzweing (2000); Rts Laney & Stobie (1995) / Rgfg Gieren et al. (1998) in the table " refer to Ru 
and * to Rgfg; Rrbmr Ripepi et al. (1997); Rbm Barnes & Moffet (1987). 



affect the radius determination. Therefore, we decided to ex- 
clude from our relations the first overtone pulsators and the 
stars flagged as "B" and "O"^ in Table |3 whereas we left in 
our sample the Cepheids flagged as "B:", "b" and "V" (i.e. 



uncertain spectroscopic binaries and separated visual binaries 
respectively, see Szabados 2003). As a result of this selection 
procedure, we are left with 20 and 16 Cepheids in the cases 



' We tried to determine the radius for the few Cepheids with known locity curve, but the error produced in this procedure remained still 
orbits, by disentangling orbital motions and pulsation in the radial ve- too high to secure good results. 



10 



A. Ruoppo et al.: Improvement of the CORS method for Cepheids radii determination based on Stromgren photometry 




Fig. 12. Comparison of radii obtained in this work {RwnhhB^ 
RnoAB) with those derived by other authors (Arellano Ferro & 
Rosenzweing (2000) RAFR;Ripe:pi et al. (1997) RRBMR.Moffett 
& Barnes (1987) 7?bm). 

without and with AB respectively. We therefore calculated new 
Period-Radius relations with the following results (see Fig.ll4>: 

log = (1.18 + 0.08) + (0.69 + 0.09) log P r.m.i = 0.08 (14) 

log/; = (1.19 + 0.09) + (0.74 ±0.11)logP r.m..s = 0.07 (15) 

A comparison between Eq. [2] and Eq. [H] ^\ shows 
that the net result of our selection criterion mainly consist in 
increasing the slope of the PR relation, going in the direction 
to improve the agreement with the literature, and, in particular 
with the NIR results by Laney & Stobie (1995) and Gieren et 
al. (1998) (see table|6lfor a comparison between present results 
and other relations in the literature). Yet, we have to note that 
the errors of the derived coefficients are rather large. This is 
caused by: 1) the small number of Cepheids left after the se- 
lection process (In particular the number of binary stars in our 
sample is very large) 2) the lack of long period Cepheids in 
the case with AB, as a result of problems with synthetic model 
atmosphere grids (see Sect. 3.1). We remark, however, that in- 
cluding or not the AB term in the CORS determination makes 
some difference, in the sense that, apart from a few Cepheids 
in the low period range, the scatter in the PR relation is slightly 
reduced (see Fig.ll4> when the AB term is included. This result 
seems to suggest that the modification to the CORS method 
presented in this paper is well suited also when applied to real 
stars (not merely to synthetic light curves). 

Finally, we notice that our results concerning the radii and 
the PR relations could be, in principle, combined with an ef- 
fective temperature calibration to derive the intrinsic stellar lu- 
minosity, through the Stefan-Boltzmann law, and in turn dis- 
tance to studied Cepheids. This possibility will be investigated 



Fig. 13. Top: Period-Radius relation obtained in the case with- 
out AB; Bottom: as before, but with AB. Open circle: Binary 
stars; Triangles: first overtone pulsators; Crosses: stars with 
loop in the [mi], [ci] plane not completely within the theoreti- 
cal grids, see par.3 

in a forthcoming paper. In the following section we will apply 
an alternative method to derive the luminosity and distance by 
means of the comparison of empirical light and radial velocity 



2.5 - 




0.5 1 1.5 

logP 

Fig. 14. Top: The same of figurefOlbut excluding first overtone 
and binary stars. 
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Table 6. Comparison between the PR (log R - a log P + b) coefficients obtained in this paper and the ones based on selected 
works in literature. In particular, from left to right, we report the slope, the zero point, the source and the method adopted for 
deriving the PR. 



a 


b 


Source 


Method 


n -L. n m 


1 1 1 -I- n m 




Surf. Brightness 


0.751 ± 0.026 


1.070 ±0.008 


LS 


Surf. Brightness 


0.750 ± 0.024 


1.075 ± 0.007 


GFG 


Surf. Brightness 


0.606 ± 0.037 


1.263 ± 0.033 


RBMR 


CORS 


0.655 ± 0.006 


1.188 ± 0.008 


BCM 


Theory 


0.69 ± 0.09 


1.18 ±0.08 


This work 


new CORS (without AB) 


0.74 ±0.11 


1.19 ±0.09 


This Work 


new CORS (with AS) 
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Arellano Ferro & Rosenzweig (2000, AFR);Bono et al. (1998, BCM). 



curves with the predictions of the nonlinear convective pulsa- 
tion models discussed above. 

6.1. Theoretical Fit of the light and radial velocity 
curves of the Cepheid Y Oph 

Before going to the conclusions of this paper, we note that it 
has recently been suggested that nonlinear pulsation models 
provide a direct tool to evaluate the intrinsic stellar properties 
of pulsating stars through the comparison of observed and 
predicted variations of relevant parameters along a pulsation 
cycle (see Wood, Arnold & Sebo 1997, Bono, Castellani & 
Marconi 2000, 2002). This kind of analysis also provides 
an additional test for our radius determination technique by 
means of the comparison with the radius of the model which 
is able to simultaneously reproduce the period, the amplitude 
and the morphology of light and radial velocity curves. We 
plan to apply this method to a sample of Galactic Cepheids 
with accurate photometric data and available radial velocity 
information. In this paper we present a first application to 
the Cepheid Y Oph. The observed properties of this star are 
summarized in Table 




PHASE 



Fig. 15. Top panel: empirical light curve (dots) for the star Y Oph, 
compared with the best fit model (solid lines, see text for details). 
Bottom panel: the same but for the radial velocity curve. 
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7 Mo 


0.28 


0.02 


4720 K 


97.77 


3.64 



Table 7. Top: Observed properties of Y Oph; Bottom: Physical pa- 
rameters of the best fit model 



Starting from the observed radial velocity and light curve, 
we try to reproduce their morphology and amplitude, by com- 
puting pulsation models along isoperiodic (with period equal to 
the observed one) sequences with varying pulsation mass and 
effective temperature and assuming, for each mass, a canoni- 
cal mass-luminosity relation. The best fit model resulting from 
these computations is shown in Fig.^](solid line). The corre- 



sponding stellar parameters are reported in Table Q and show 
that the agreement with the radius determined in this paper is 
good. 

The other important information provided by the fit is the 
stellar absolute magnitude My - -3.996 mag, which allows 
us to estimate a distance of about 423 pc, consistent with the 
independent evaluation by Gieren, Barnes & Moffett (1993), 
thus supporting the predictive capability of pulsation models. 

7. Final remarks 

We have presented a modified version of the CORS method 
based on a new calibration of the Surface Brightness function 
in the Stromgren photometric system. In particular we have 
been able to derive a calibration of 5 y as a function of the 
Stromgren reddening free indexes [mi] and [ci] by adopt- 
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ing grids of theoretical colors. This procedure revealed the 
unexpected occurrence that the quoted theoretical grids are 
not able to fully include the location of actual long period 
(P > 12-13 days) Cepheids loops in the [mi],[ci] plane. This 
problem could only be overcome by adopting next generation, 
hopefully improved, model atmosphere grids. Nevertheless, 
the modified CORS method presented here has been tested 
by means of synthetic light and radial velocity curves derived 
from nonlinear pulsation models and simulations have been 
performed to take into account the quality of real observed 
curves as well as possible time shifts between photometric 
and radial velocity data. The results of such tests can be 
summarized as follows: 

a) the present method appears capable to derive the radius of a 
Cepheids with an average precision of around 3.5% and 5% in 
the case of good and fair observational data respectively (see 
Section 4.3); 

b) the error associated with a phase shift between photometric 
and radial velocity data could be as large as ~ 17% in the 
extreme case of ~ 0.05. In a more common case of 
6.(p ~ 0.01 the relative error is still rather large, namely ~ 5%. 
If we add to the error budget the uncertainty on the value of 
the projection factor k {~ ±2% for M ± 0.03), neglecting 
other possible contributions (see Bono, Caputo & Stellingwerf 
1994 and references therein for a detailed discussion), we 
can conclude that estimating the radius of a Cepheid by 
means of the Baade-Wesselink method (at least in the present 
formulation) with an accuracy better than 10% is extremely 
difficult to be achieved. 

The method has been then applied to a sample of Galactic 
Cepheids with Stromgren photometry and radial velocity data 
to derive the radii and a new PR relation. As a result we ob- 
tained the following Period-Radius relation: log/? - (1.19 + 
0.09) -I- (0.74 ± 0.11)logP (r.m.s=0.07). This relation, even if 
not very accurate (mainly due to the unexpected presence of a 
large fraction of binary stars in our Cepheid sample), neverthe- 
less is in satisfactory agreement with previous findings in the 
Uterature. 
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Appendix A: 

The 4"" degree polynomial fit to effective temperature, effective 
gravity, and bolometric corrections mentioned in Sec. 3. 1.1 are 
the following: 

log Teff = ao + a\- [mi] + ■ [rriif + ■ [ci] (A.l) 
+ fl4 ■ [mi] ■ [ci] + fl5 ■ \m\f ■ [ci] +a(, ■ \c\f 
+ ai ■ [m\f ■ [c\f 

log geff = bo + bi- [mi] + b2 ■ \mif + bj, ■ [ci] (A.2) 
+ bi ■ [mi] • [ci] + bi ■ [cif + be ■ [mxf ■ [cif 



BC = Co- log Teff + Cl • log T^ff (A.3) 

+ C2 ■ log geff + C3 • log Teff ■ log geff 
+ C4 ■ log T^fj: ■ log geff + C5 • log gljf 
+ C6 • log Teff ■ log glff + Cj ■ log T^ff 

the coefficients a,, bj, c, of the previous relations are listened 
in Table A.l and the other symbols have their usual meaning. 
Note the r.m.s of the previous relations are 0.0018 dex. , 0.1 
dex. and 0.003 mag. respectively. 



Table A.l. Coefficients for the polynomial fits described in the 
Appendix. 
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